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We have studied the sensitivity of the ATLAS detector for supersymmetric neutralino signals in the Zj + 
final state in a GMSB model in which the Higgsino-like neutralino is NLSP. This study considers the reaction of 
PP ~* XiXi ~* ^{~^ £~^£~)'jGG, where £ = /i/e, at a center of mass energy of y^=10 TeV using fully simulated 
ATLAS Monte Carlo events both for the signal and background. Based on the GMSB Model Line E predictions, 
we expect that, for an integrated luminosity of less then 3 fb~^, ATLAS could detect the GMSB signal from 
the Higgsino-like neutralino that has a mass of 134.7 GeV in the Z7 + final state with a significance of 5cr, 
assuming 20% systematic uncertainty. 



1. Introduction 

One of the major motivations of the LHC exper- 
iments is to search for physics beyond the Standard 
Model (SM). Supersymmetry (SUSY) is a favored can- 
didate in this search. However, SUSY must be bro- 
ken so that the masses of super particles exceed those 
of their SM partners. The minimal supersymmetric 
extension of the standard model (MSSM) does not 
incorporate SUSY breaking and has 124 free param- 
eters. SUSY breaking can be achieved by introduc- 
ing a hidden sector. The mediation of the break- 
ing from the hidden sector to the observable MSSM 
sector is enabled by different mechanisms such as 
mSUGRA, gauge-mediated SUSY breaking (GMSB) 
and AMSB through gravity, gauge or anomalous in- 
teractions. Through this mediation, the number of 
free parameters is significantly reduced. 

We have investigated neutralino pair production at 
the LHC at a center of mass energy of ^/s — 10 TeV 
within the framework of a GMSB scenario [l|. A dis- 
tinctive phenomenological feature of the GMSB model 
is the presence of neutral dibosons {Z or 7) plus miss- 
ing transverse energy (^t) in the final state. 

This note reports on the sensitivity of the ATLAS 
detector at ^/s = 10 TeV for the process pp XiXi 
Z{-^ i^£~)jGG, where £ — ^/e. Xi in this process 
is a Higgsino-like neutralino, which is the next-to- 
lightest supersymmetric particle (NLSP) in the model 
considered here . It decays to a neutral vector boson 
{Z or 7) and a gravitino (G), the hghtest supersym- 
metric particle (LSP), which escapes detection result- 
ing in large ■ 

1.1. Theoretical overview 

In GMSB models, SUSY, which is broken at a 
TeV scale in a hidden sector, is propagated down 
to the MSSM observable sector via new chiral su- 
permultiplets, called messengers, that couple indi- 
rectly to the MSSM particles through the ordinary 



SU{3)c X SU{2)l X U{1)y gauge boson and gaug- 
ino interactions. The main advantage of GMSB mod- 
els is the automatic creation of identical soft SUSY 
breaking masses for scalars with the same gauge quan- 
tum numbers but different flavors. Therefore, there 
are no problems with flavor changing neutral currents 
(FCNC) or CP violation constraints. 

Within the minimal version of the GMSB frame- 
work (mCMSB), the couplings, branching ratios, de- 
cay widths and sparticles masses are determined by 
the parameters: 

A, M, TVs, ian{l3),sign{^i),Cgrav (1) 

A is the SUSY breaking energy scale, A = F/M, 
where is a vacuum expectation value of an auxiliary 
field that determines the magnitude of supersymme- 
try breaking in the vacuum state. M is the size of the 
messenger mass scale, M > A. For electroweak scale 
superpartners, A is ~ 100 T€\f /^/N^. is the num- 
ber of generations of messenger fields, tan(/3) is the 
ratio of the MSSM Higgs vacuum expectation values 
{{Hu)/{Hj)), sign{^) is the sign of the Higgs sector 
mixing parameter and Cgrav is the ratio of the mes- 
senger sector SUSY breaking order parameter to the 
intrinsic SUSY breaking order parameter, which con- 
trols the coupling to the gravitino. 

In GMSB, G is the LSP (usually, M(G) << 1 GeV). 
The G has a derivative coupling to each particle and 
its superpartner with an interaction strength that is 
inversely proportional to V^. Because of this cou- 
pling, the next-to-lightest superpartner is unstable 
and decays to its lighter partner through G emission. 
As a result, the nature of the NLSP defines the phe- 
nomenology of the GMSB model. Depending on the 
region in parameter space, either the neutralino (x^) 
or the f arise as the NLSP. The NLSP decay length, 
which depends on ^/F, can be divided into three 
ranges: (1) the NLSP decays promptly, (2) the NLSP 
decays inside the detector away from the collision 
point {VF < 10^ GeV), and (3) the NLSP decays 
outside the detector (a/F > 10^ GeV). We consider 
only the case where the NLSP is the Xi, which decays 
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promptly, corresponding to Cgrav ~ 1- 

x1 is a mixture of gaugino {B, W'^) and Higgsino 

{H^, H^) eigenstates, and therefore decays to a 
7, Z, or Higgs (h). If Xi is gaugino-like, it decays 
mostly to jG, leading to the 77 + signature (with 
R-parity conservation, supersymmetric particles are 
produced in pairs). If the Xi is Higgsino-like, it de- 
cays to hG. In addition, since the longitudinal compo- 
nent of the Z boson mixes with the Goldstone mode 
of the Higgs field, a Higgsino-like neutralino also de- 
cays to ZG. Because of a strong phase space sup- 
pression of the h and Z final states, decay to a 7 can 
also be important for Higgsino-like neutralinos, which 
are not very much heavier than the Z boson. Conse- 
quently, a pair of Higgsino-like neutralinos produced 
in a collider can give rise to the diboson final states 
{jj,hh,h'y,hZ,Z-/,ZZ) + 3, H- We consider 
the possibility of a Higgsino-like xf NLSP, where one 
neutralino decays to -fG and the other decays to ZG, 
leading to the Z-f + signature. 

1 .2. Z7 + signature 

To study the Z7 + signature in the GMSB 
framework, we generated a sample of pp — > XiXi ~^ 
Z{£'^£~)^GG events, where £ = fi/e, using a 
(non-minimal GMSB) Higgsino-like neutralino model 
(model line E), described in This model is char- 
acterized by the following parameters: 

M/A = 3, iVs = 2, tan(/3) = 3, = 0.75Mi, (2) 

where Mi is the gaugino mass and yields 225.4 GeV. 
Unlike in the mGMSB, where the absolute value of /i 
is set by the electroweak symmetry breaking (EWSB) 
radiative conditions as shown in equation [3], ^ is set 
to O.75M1 in this approach, which gives rise to a large 
Higgsino fraction in the neutralino. 

^ + tan^(/3)-l 

Figure [H taken from 0, shows the production 
cross-sections for sparticles(left plot) and branching 
fractions of Xi (right plot) in model line E, as a func- 
tion of the most influential parameter, A. The dom- 
inant source of sparticles in this model is produc- 
tion of neutralino pairs, chargino pairs {xf, i = 1,2), 
and chargino-neutralino pairs. Charginos and heav- 
ier neutralinos decay to Xi through cascade decays. 
Gluinos and squarks are relatively heavy; therefore, 
their production becomes suppressed when A is in- 
creased The branching ratio Br(xi ZG) be- 
comes dominant for A > 80 TeV (M(x?)~134 GeV); 
we thus choose this value of A for the current study. 
This choice gives a production cross-section of the de- 
cay XiXi ~^ Z{1^ £^)^GG , where t — /i/e, of about 



56.2 fb. The cross-section is calculated as follows: 
2 X a{pp ^ x?X?) X Br(x? ^ ZG) x Br(x; ^ 7G) 
X Br(Z ^ l+fr), where (T{pp x?X?) = 1880 fb, 
Br(xi 7G) = 0.55, Br(x? ^ ZG) = 0.4, and 
Br(Z l+i-) = 0.067 [|. The SUSY particle mass 
spectrum for model line E with A = 80 TeV is listed 
in Table H 



Table I SUSY particle mass spectrum for the GMSB 
model characterized by the following parameters: A = 
80 TeV, M = 240 TeV, V5 = 2, tan(/3) = 3, = 
169 GeV, M{t) = 172.5 GeV. 



Names 


Mass Eigenstates 


M (GeV) 




Ul 


1272.9 




Ur 


1220.6 




dL 


1272.9 




dR 


1215.6 


squarks 


Sl 
Sr 


1272.9 
1215.6 




Cl 


1272.9 




cr 


1220.6 




ii 


1119.1 






1239.2 




bi 


1216.9 




b2 


1228.4 




Cl 


408.9 






196.1 


sleptons 


p.L 


408.9 




fiR 


196.1 




fl 


198.9 






407.9 




X? 


134.7 


neutralinos 


5& 

x§ 


172.3 
241.3 




xl 


443.3 


charginos 


xf 
xf 


152.1 
442.9 


gluino 


9 


1221.1 


gravitino 


G 


4.6x10"^ 






97.4 


Higgs bosons 


H° 

A° 


746.5 
739.3 






748.5 



2. Signal MC production and decay 

The mass spectrum and branching fractions were 
calculated using ISASUGRA 7.75 and ISASUSY 
7.75 @. ISASUGRA takes input parameters from 
a particular SUSY model. We used the follow- 
ing mGMSB parameters: A = 80 TeV, M = 240 TeV, 
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Figure 1: Cross-sections for various sparticle production processes in model line E (left). The branching fractions for 
various decays of the Xi NLSP (right) in model line E versus the parameter A. The scale on top shows the Xi mass. 



= 2, tan(/3) = 3, sgn{^) — 1, Cgrav = 1 to pro- 
duce the corresponding weak scale MSSM parameter 
set. One of these parameters, /z, was then changed ac- 
cordingly to the model we used: fj, = 0.75Mi. These 
parameters were then taken by ISASUSY to evalu- 
ate the sparticle masses, decay rates, and branching 
fractions that serve as input for the program which 
generates the full SUSY event. About 50k events of 
XiXi ~^ Z{£'^£~)jGG, where i = ji/e, were gener- 
ated using HERWIG/ JIMMY 0,9, which generates 
the sparticle cascade decays, parton showers, hadro- 
nisations, and underlying events. The full ATLAS 
detector simulation was done using GEANT4 % All 
signal and background MC events were reconstructed 
using standard ATLAS reconstruction software. The 
geometrical acceptance cuts (|??(^)| < 2.5) reduced the 
production cross-section of x\Xi Z{£~^£~)jGG to 
about 41 fb. 



3. Background 

The major background processes are summarized 
in Table |TT1 We consider three types of background: 
(1) real Zj pairs from SM Z{^ £~^£~) + X processes 
which include initial state radiation (ISR) of photons 
from the colliding quarks and final state photon radi- 
ation (FSR) from the Z decay leptons; (2) events with 
one real Z and jets which are misidentified as 7's; and 
(3) top pairs which decay to leptons: tt £ + X. 

4. Event selection 

To identify muons, electrons, and photons we used 
the following criteria: 



Table II Background NLO production cross-section after 
event selection {(Jnlo), fuUy simulated number of MC 
events (Nmc), and the luminosity which corresponds to 
the MC events. The last column indicates the MC gener- 
ators used to produce the MC events. 



Process 


onlo 
(Pb) 


Nmc 


L 

(fb-^) 


Generator 




7.1 


22489 


3 


PYTHIA 


tt^e + x 


202.9 


1860622 


9 


MC@NLO 




1317.8 


4909710 


3.7 


PYTHIA 


Z{^ e+e-)+X 


1317.8 


4766732 


3.6 


PYTHIA 


W+ Z Ivll + X 


264.73 


20000 


75.5 


MC@NLO 


W'Z ^ £vU + X 


155.96 


15457 


99 


MCONLO 


W+Z qq££ + X 


828.52 


5000 


6 


MC@NLO 


W'Z qq££ + X 


488.1 


5000 


10 


MC@NLO 


W+Z Tu££ + X 


132.37 


19675 


149 


MC@NLO 


W-Z Tu££ + X 


77.98 


19719 


253 


MC@NLO 



• Muons are identified with an algorithm which 
associates a track found in the muon spectrom- 
eter, after corrections for energy loss in the 
calorimeter, with the corresponding inner detec- 
tor track. This association is performed by sta- 
tistically combining the two fitted sets of track 
parameters using their corresponding error ma- 
trices. The rapidity acceptance of the combined 
muon track is limited by the inner detector to 
I77I <2.5. 

• Electrons are identified using a multivariate 
Boosted Decision Trees (BDT) classification al- 
gorithm. BDT input variables are selected from 
measurements in the hadron calorimeter, the 
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electromagnetic calorimeter and the inner de- 
tector. In the hadron calorimeter the principal 
variable is the fraction of hadronic Ex leakage 
of the total cluster Et- The first two longitudi- 
nal samplings of the electromagnetic calorimeter 
provide sets of correlated variables, including, 
for example, the fraction of energy deposited in 
each sampling. Track parameters and the num- 
ber of hits in each of the inner detector layers 
provide a third set of variables; and finally, a set 
of inner detector track and EM cluster matching 
parameters are used. There are approximately 
20 input variables incorporated by the BDT al- 
gorithm. 

• Photons are selected using the 'tight' category 
of photon identification cuts described in 

To select events with ^+^~) and 7, the fol- 

lowing preselection cuts are applied. To select a 
good lepton candidate, the pt{^) is required to be 
greater than 6 GeV with \rj\ < 2.5. A cut of 30 
GeV is applied on the leading lepton. In the pre- 
selected lepton sample, we search for di-lepton can- 
didates with 70 GeV < M {£+(-) < 100 GeV. To se- 
lect 7 candidates, the pril) is required to be greater 
than 20 GeV. To reject 7 produced through lepton 
bremsstrahlung radiation, the minimal value of AR, 
where AR = y/Aif + AcI^, between selected leptons 
and 7 is required to be greater than 0.1. Table UTTl 
shows the efficiencies and effective remaining cross- 
sections after preselection cuts are applied to the sig- 
nal and background processes. Figure [2] shows the in- 
variant mass distribution of £'^£~' after preselection. 

4.1. + selection 

After preselection, significant background contam- 
ination remains from SM ^(— > £~^£~)j and inclusive 
Z{-^ £~^£~) + X processes, which have large produc- 
tion cross-sections. The former has a topology sim- 
ilar to the signal, where a real Z and 7 tend to be 
back-to-back in the transverse plane. However, the 
reconstruction efficiency of this sample is one order of 
magnitude lower than the signal efficiency due to the 
Pt{i) cut, applied in the preselection. For the lat- 
ter, the source of the 7's is mostly misidentified jets; 
therefore, the reconstruction efficiency is low. For fur- 
ther rejection of background, various distributions of 
kinematic variables of the background and signal MC 
samples were compared. 

The lepton pair is required to have an invari- 
ant mass close to the Z mass, specifically \Mii — 
91.2 GeV| < 13 GeV. This requirement is equivalent 
to a 5(7 width, and helps to reduce background where 
the lepton pair does not come directly from a Z . 

The minimal value of AR between leptons and 7 
is required to be greater than 0.4. This cut rejects 



Table III Results after preselection for the signal and 
background MC samples. The number of preselected 
events (Nps), the efficiency of the preselection (tpa) and 
finally the cross-section after preselection {(Jps) are given. 





Process 


Nps 


€ps (%) 


aps (fb) 




X°iXi ^ fi+n' 




11574 


49 


10 




Z{^ £+£-)j 




1845 


8.2 


581 




tt^e + x 




127 


0.007 


14 




Z{^H+^i')+X 




2440 


0.05 


655 




Z{^ e+e~)+X 
W+Z lv££ + X 






168 


0.8 


2 




WZ ^ £m££ + X 




162 


1 


1.6 




W+Z qq££ + X 




8 


0.16 


1.3 




W'Z qq££ + X 




8 


0.16 


0.8 




W+Z Tu££ + X 




55 


0.28 


0.4 




W'Z ^ Tv££ + X 




55 


0.28 


0.2 




Xlx^i - Z{~. e+e- 


)lGG 


6985 


30 


6 




Z{-^ £+£-yf 




1187 


5.3 


374 




tt^e + x 




25 


0.001 


2.7 




Z{^ ^i+fi-)+X 











2e 


Z{^ e+e-)+X 
W+Z ^ lv££ + X 
W'Z -*£v££ + X 
W+Z qq££ + X 




1263 
116 
88 
9 


0.026 
0.6 
0.6 
0.2 


349 
1.5 
0.9 
1.5 




W'Z qq££ + X 




7 


0.14 


0.7 




W+Z Tu££ + X 




38 


0.2 


0.25 




W~Z -> Tiy££ + X 




42 


0.2 


0.17 



mostly tt and ^(— > £+£-) + X events where the real 
7, produced through lepton bremsstrahlung radiation, 
is in the same direction as leptons. 

Signal events are characterized by the presence of 
relatively large ^t- To reject the SM diboson pro- 
cesses that do not have sources of significant ^t, we 
required ^Jt > 40 GeV. Figure [3] shows the f^r dis- 
tributions of the signal and background MC events 
normalized to the cross-section after preselection. 

The transverse like mass quantity (Mrp) 10] of two 
leptons, 7's, and ^t, defined as 



(4) 

is highly correlated with ^t- Unlike the ^t, Mj, helps 
reject tt events. We required Mj, to be greater than 
210 GeV. 

The Zj system from the decays Xi iG and x? ^ 
ZG is strongly boosted by the neutralinos and has a 
high invariant mass and px relative to the Z7 from 
the decay of Z{-^ £+£') + X . The M(Z7) is thus 
required to be greater than 110 GeV. 

The photons from decays of Xj* are strongly boosted 
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. ATLAS Work in Progress 



ZHnH) + X 




o 0.2- 



0.1 



ATLAS Work in Progress 



X°Z° ^ Z(ee)YGG 




100 



M [GeV] 



M [GeV] 



Figure 2: Distribution of M(/i^/i ) (left) and M(e'''e ) (right) after preselection. The plots are normalized to unit area. 




Figure 3: Distributions of after Z{^'^ ^ )7 (left) and Z(e"'"e )7 (right) preselections. The plots are normalized to 
cross-section. 



and have a high relative to the fake photons and 
ISR photons, which populate the low px region. The 
Pt of the 7 is required to be greater than 30 GeV. 
The Pt of muons and electrons from the decay x\ —> 
Z{i'^£~)G is high relative to the low pr leptons which 
decay from b hadrons. Therefore the cut of 10 GeV 
applied to the minimum px leptons helps to reject the 
tt background events. 



In addition, ti and Z{-^ t^l^ ) + X events can be 
effectively reduced by track isolation conditions on the 
leptons and photons. The isolation condition consists 
of requiring less than four tracks in a cone of Ai? = 
0.3 around the reconstructed lepton and 7. A cut of 
Ado < 0.1 mm (Ado < 0.15 mm) for muons (elec- 
trons) and Azo < 1 mm, where Ado(Azo) is the dif- 
ference between the transverse(longitudinal) impact 
parameters of the tracks, helps to reject tracks from h 
hadron decays which are expected to have, on average, 
large transverse and longitudinal impact parameters. 



5. Results 

Table IIVI summarizes the results of the analysis 
for an integrated luminosity of 1 fb~^. Since the 
Z{-^ t^l~) + X sample includes ISR photons, the 
Z{-^ t^l^)^ sample was not included in the total 
background estimation to avoid the double counting 
of Z7 events. After selection, four events remain in the 
Z(— > [L^ pr)+X sample; two of these events (0.5±0.4 
events for 1 fb"^) include ISR photons. This number 
is consistent with the Z(^-^ t^t')^ sample, which has 
one event (0.3±0.3 events for 1 fb"^) that survives the 
selection criteria. The two remaining events surviving 
selection cuts from the Z{^ [i^ pT) + X sample in- 
clude fake photons. The number of expected events 
from the 2^(— > A*^/^" ) + ^ sample due to jets misiden- 
tified as photons is thus expected to be 0.5 ± 0.4. The 
dominant contribution of the background comes from 
Z{-^ ^i+fi-)+X (24%) and Z{-^ i+£-)-f (24%) sam- 
ples after muon selection and from tt (40%) and in- 
clusive Z (30%) samples after electron selection. 
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Table IV Results of cut-based analysis. This table in- 
cludes the number of expected events for 1 fb~^. The 
errors shown are statistical or 90% CL in the case of in- 
sufficient statistics. 



2/i 



2e 



Process 



XiXi ^ 
tt-*£ + X 

ivll + X 



W+Z 
W~Z 
W+Z 

w-z 

W+Z 



■ £u££ + X 
qqll + X 

■ qqU + X 
TvU + X 



W'Z TvU + X 
Total bkg {Ni) 



z{-^ rr)'y 
tt^i + x 

Z{^ e+e-)+X 
W+Z ^ IvU + X 
W'Z ^ £uU + X 
W+Z -> qqll + X 
W'Z -* qqii + X 
W+Z ^ TuU + X 
W'Z TuU + X 
Total bkg {Nb) 



N of expected events 
(for 1 fb"^) 



5.4 ± 0.06 
0.3 ± 0.3 
0.3 ± 0.2 
1.1 ± 0.5 

0.28 ± 0.06 
0.26 ± 0.05 

< 0.4 

< 0.2 
0.07 ± 0.02 
0.04 ± 0.01 

2.05 ± 0.6 



3.38 ± 0.05 

< 0.7 
0.4 ± 0.2 
0.3 ± 0.3 

0.17 ± 0.05 
0.06 ± 0.02 

< 0.4 

< 0.2 
0.06 ± 0.02 
0.02 ± 0.01 
1.0 ± 0.35 



for the 2e sample, assuming 20% systematic micer- 
tainty (Table W} . The overall detection significance 
of the XiXi ~^ Z(t+£~)^GG decay, where £ = ^/e, is 
expected to be 3.2cr assuming 20% systematic uncer- 
tainty [l^l for an integrated luminosity of 1 fb~^. For 
an integrated luminosity of 3 fb~^ the significance is 
expected to be 4.2cr for muons, 3.6f7 for electrons, and 
5.6cr for the combined samples assuming 20% system- 
atic uncertainty. 



6. Summary 

This note presents a study of a GMSB search at 
ATLAS with a Z7 -I- final state where the Z is 
reconstructed using ee or /i/i based on fully simulated 
MC events. The study has been performed within 
the GMSB model hue E with A=80 TeV, using 5_0k 
fully simulated signal events of XiXi — ^ Z{£+£~)jGG 
that are characterized by two high px leptons plus 
an energetic 7 with fix- The results show that a 
significance of 5.6a, assuming 20% systematic uncer- 
tainty, can be attained with an integrated luminosity 
of 3 fb""'^ in the combination of the Z{e+e~)^ and 
Z{p+ final states using a cut-based method af- 
ter taking into account the background contributions 
from SM ti, Z{~^ £+£-)+ X , and WZ. This suggests 
the feasibility of a search for GMSB at ATLAS with 
3 fb~^ of integrated luminosity given a center of mass 
energy = 10 TeV. 
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z{^^+^i-)J + 


^(e+e-)7 +-^T 


combined 
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0.04 
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Significance 
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l.Scr 


3.2(7 
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